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, Using perturbative methods, we analyse a nonlinear generalisation 

°f Schrodinger's equation that had previously been obtained through 
information-theoretic arguments. We obtain analytical expressions for 
the leading correction, in terms of the nonlinearity scale, to the en- 
ergy eigenvalues of the linear Schrodinger equation in the presence 
£3 ■ of an external potential and observe some generic features. In one 

space dimension these are: (i) For nodeless ground states, the energy 
q-( shifts are subleading in the nonlinearity parameter compared to the 

shifts for the excited states, (ii) the shifts for the excited states are 
due predominantly to contribution from the nodes of the unperturbed 
wavefunctions and (iii) the energy shifts for excited states are posi- 
tive for small values of a regulating parameter and negative at large 
values, vanishing at a universal critical value that is not manifest in 
the equation. Some of these features hold true for higher dimensional 
problems. We also study two exactly solved nonlinear Schrodinger 
equations so as to contrast our observations. Finally, we comment on 
the possible significance of our results if the nonlinearity is physically 
realised. 
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1 Introduction 



Various nonlinear extensions of Schrodinger's equation have been proposed 
[lj over the years as possible generalisations of the linear evolution of the 
original theory. Although several low-energy experiments have placed very 
small upper bounds [2] on the proposed extensions, there is still the possibility 
that quantum mechanics might have to be modified at high energies or short 
distances [3] where the structure of spacetime is expected to be different [I] . 

However in this paper we remain within the non-relativistic realm so as 
to explore in more detail the properties of one particular nonlinear extension 
that was motivated in [5] by maximum uncertainty arguments [HIE] similar to 
those used in statistical mechanics [8]. In higher than one space dimension the 
equation of Ref. [5] was not rotationally invariant, motivating a suggestive link 
between spacetime symmetries and quantum linearity. Some implications of 
such a connection for phenomenology were discussed heuristically in [5l [3] . 

In Ref. [5] some exact nonperturbative solutions of the abovementioned 
equation were obtained, displaying intriguing and novel features that are 
probably related to the unusual structure of that equation. Indeed, it was 
hinted in [9] that the equation might also be interesting as an effective equa- 
tion in other domains of physics, such as nonlinear optics [10], rather than 
its original intention in Ref. [9]. 

Here we investigate how the nonlinearity of that equation perturbs the 
energy spectrum of the usual linear Schrodinger equation. Since simple esti- 
mates already indicate that the size of the nonlinearity scale must be tiny for 
it to be consistent with phenomenology [5], we shall use standard first order 
perturbation theory for our study. Our primary aim here is not to confront 
empirical data but to uncover further properties of the nonlinear equation. 
As we shall see, even at the perturbative level the equation of Ref. [5] has a 
rather surprising character. In particular, we find a universal critical point 
of the theory that is not at all obvious from the equations of motion. 

In the next section we outline our perturbative scheme and then illustrate 
it with some numerical results in Section (3). A general analytical investi- 
gation is next conducted in Section(4) to extract and explain the features 
mentioned in the abstract. In Section (5) we contrast the perturbative prop- 
erties of the nonlinear equation with those of two other exactly solved non- 
linear Schrodinger equations. The concluding section summarises the main 
lessons and discusses some implications. The appendices contain additional 
derivations. 
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2 Perturbative framework 



Let us focus first on the nonlinear equation for a single particle in one space 
dimension that was derived in Ref.[5], 

with p(x,t) = $?*(x, t)if?(x, t) the conserved probability density and 

F(p) = Q n i-Q, (2) 

where 



Qn 



h 
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rj)p + r/p + (1 — r])p + r]p + (l — rj)p_ + r)p 

(3) 

is a regularised nonlinear "quantum potential" . The parameter 7] takes values 
< r] < 1, its crucial role being to regulate potential singularities where p(x) 
vanishes. We have used the notation 

P±( x ) =p(x±r]L). (4) 

Note that if \l/ is any solution of the equation, then so is A\l/ for an arbitrary 
constant A, so we may re-normalise states freely. The nonlinearity is charac- 
terised by the length scale L, in terms of which one may perform a formal 
expansion of ([3]), 

2m ^Jp dx 2 ' 

with a remainder of O(L). 

Let * = e- iEt/h (f){x) be the energy eigenstates of the usual linear Schrodinger 
equation for a given external potential V(x). Assuming that the spectrum 
deforms continuously as the nonlinearity F is turned on, then to leading 
order the corrected energies are given by first-order perturbation theory, 

E exact = E + 5E, (6) 

rod 

5E = dx 4>*F(<j>)<j> . (7) 



Note that F is evaluated using the unperturbed wavefunctions and so from 
now on p will refer to <f>*(x)<f>(x). 
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We may trust first-order perturbation theory when the nonlinearity is 
small. The relevant dimensionless expansion parameter is L/a where a is 
a typical scale in the linear theory, such as the deBroglie wavelength. As 
equations fl3|4|) indicate, the expression ([7j) is actually a complicated function 
of L/a from which the leading behaviour must be extracted. We will discard 
subleading terms from (j7j) as they will be of the same order as second-order 
perturbation theory corrections, which we do not study here. 

For the problems we will study, the unperturbed wavefunctions <f>(x) are 
parity eigenstates so that p(x) = p(—x). Changing variables x — > — x in (J7J) 
and using the parity invariance of p shows that 

6E(L) = 5E(-L) . (8) 

That is, when we allow the parameter L to take negative values, then al- 
though equation flTJ is not invariant under L — > —L, yet the first-order 
energy shifts are. Therefore if 5E(L) were an analytic function of L, one 
would have concluded that 

8E(L) ~0(L 2 ) + 0(L 4 ) + ... (9) 

as the 0(L°) term vanishes by construction, see (jSJ). In reality however, 
5E(L) is generically non-analytic! To see this, consider the naive series ex- 
pansion of the integrand in (jTJ). It results in the formal expression 

6E(L) oc LV / [6(2- 3r?)V) 4 - 12(3- 877 + 6 V 2 )p(p') 2 p" 

Joo p A 

+ Ap 2 p'p"' + p 2 {3(p") 2 - 2pp"")\ (10) 

which is ill-defined because of the singularities that occur where p{x) van- 
ishes, that is where the unperturbed wavefunction has nodes. Thus one may 
conclude SE(L) ~ 0(L 2 ) only for nodeless states, which are typically only 
the ground states of a system. 

Since excited states of the unperturbed theory have nodes, we cannot use 
f fTUj) for them. In the next section we perform a numerical investigation of 
expression ([7j) and then return in Section (4) to a more general analytical 
investigation that explains the various observed results, such as 5E(L) ~ 
0(|L|) for states with nodes. 
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3 Numerical Investigation 



The purpose of the numerical study is two-fold. Firstly it helps us uncover 
some interesting features of the complicated nonlinearity and so guides us 
in the later, more general, analytical investigations. Secondly, it will provide 
us with important checks on the analytical derivations of Section (4) and in 
particular answer the question of how small, numerically, the perturbative 
parameter L/a has to be in the analytical expressions. 

As convergence near the end-points rj = 0, 1 is slow, we integrate ([7j) 
numerically at the symmetric point rj = 1/2, defering a discussion of other 
rj values to Sect. (4). Although phenomenologically one expects L/a to be 
tiny [5], we study much larger values ~ 10~ 3 for computational efficiency 
However we do demand SE/E ~ 10~ 2 or smaller so as to be safely in the 
perturbative regime. For each V(x), we obtain the leading dependence of 5E 
on L/a and the principal quantum number. The numerical results are then 
parametrised using a best fit to simple analytical power law expressions. 

The numerical work was performed with Mathematica [UJ and the quoted 
numbers are accurate to about the last digit. 

In the numerical work we have set L — 1 to define the reference units. 
Thus l/a factors quoted below actually correspond to the dimensionless 
quantity L/a. We have checked that the numerical results are invariant 
under L — ► — L as required by the parity invariance argument of Sect. (2). 



3.1 Infinite Well 

The infinite well with walls at x = and x = a gives 

<Pn{x) = J-sin^^ (11) 
V a a 

and unperturbed energies 

In the presence of the nonlinearity the energies shift and are given to leading 
order by 

E n = E° n + 5E n . (13) 

It is convenient to define dimensonless quantities by dividing the above equa- 
tion by h 2 7r 2 /2ma 2 , 

E n = n 2 + 5E n . (14) 
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For various fixed values of n, 1 < n < 50, the energy shifts were evaluated 
numerically for 1000 < a < 10, 000. Fig.(l) shows a log-log plot for the n = 1 
case from which one deduces 5E = — 0.99/a. The other n values give similar 
plots, all indicating 5E oc —1/a. 

On the other hand, for various fixed a, an evaluation over the range 
5 < n < 50 shows 5E oc — n 3 . Re-inserting l U we find, averaging the best fit 
for various a values, 

5E n = -1.03^-^i + 0(L/a) 2 . (15) 

In Appendix A we will explain this result analytically. 

Notice that the correction (TT5]) grows with n and so at some large value 
of n it is no longer small compared to the unperturbed value. This simply 
means that one must then go beyond first order perturbation theory. We 
discuss the possibilities in the concluding section. 

3.2 Simple Harmonic Oscillator (SHO) 

The potential is now V(x) = kx 2 /2 giving the usual unperturbed wavefunc- 
tions P2] 

Mx) = (™ 2 )~ 1/4 "° © exp (16) 

and unperturbed energies 

^ [n+ \ )Tl {i- (17) 

The linear length scale 'a' in this problem is the deBroglie length h 1 ^ 2 / (km) 1 / 4 . 
We investigated the dimensionless energies shifts 

numerically over the range, < n < 18 and 100 < a < 1000. 

For the ground state, n = 0, we find 5E oc —1/a 2 , which is a faster 
drop than seen for the infinite well. However excited states have a similar 
behaviour in 'a' to those of the infinite well, with 5E oc —1/a for any fixed n. 
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The n— dependence for fixed a is more complicated as indicated in Fig. (2). 
In summary we find for the excited states, n > 1, 



n 1A1 \L\ 

5E n = -0.2Q l -^ + 0(L/a) 2 . (19) 

a 

It must be emphasized that the result (fT9"|) is a best fit to an assumed 
power law over the limited range investigated. However, independent ana- 
lytical estimates in Appendix B do give a similar result over the same range. 



3.3 Hydrogen Atom 



We use the standard unperturbed wavefunctions as given, for example, in 
Ref.l 



lpnlm(r, 9, (f)) 



(20) 



with p = 2r/na, and the corresponding unperturbed energies 

n 2 



El 



2Ma 2 n 2 ' 



(21) 



where M is the electron mass. 

The three dimensional version of eq.([2]) is [5] 



F(p) = Q 3 -Q, 

3 h 2 
Q 3 = E 



(22) 



In 



V 



Q 



\ AML 2 r) 4 [ (1-T])p + r)p + , 
h 2 \2didip dipdip 



+ 1 



1 - T))p 



Tjp_ 



1 - rj)p + r]p +i (1 - rj)p_i + rjp_ 



8M 



P 



(23) 



with i = 1, 2, 3 and p±i(x) = ± ?7L, x 2 , x 3 ) and so on. The Bohr radius 
defines a = h 2 /Me 2 for this problem and the dimensionless energy shifts are 



6E n = 6 -^(2Ma 2 ). 



(24) 



Note that the nonlinearity breaks rotational invariance in the above ex- 
pression (|2"2"|) which is defined in the preferred Cartesian basis as discussed 
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in Ref. [5]. Thus the wavefunctions (120]) are first converted to the Cartesian 
basis for use in (I22p but the final numerical integration was performed after 
converting back to spherical coordinates. We used the built-in Monte Carlo 
subroutine in Mathematica [11] for this case and investigated only a very 
limited range of parameter values due to the time-intensive nature of the 
three dimensional problem. 

Although the pure Coulombic hydrogen atom states have a degenerate 
spectrum, we still used the simple non-degenerate first order perturbation 
theory formula for all states as our primary objective is to observe the effects 
of the nonlinearity on energy shifts. (Note also that the energy shifts due 
to the nonlinearity are expected to be much less than other effects, such as 
relativity, that in reality lift the degeneracy of the unperturbed states.) 

Consider first the zero angular momentum, I = states. For the ground 
state, u = 1 we found 6E oc — 1/a 2 while for the n = 2, 3 excited states 
we have 5E oc —1/a. This dependence on a is similar to that of the SHO. 
The dependence of the energy shifts on the principal quantum number how- 
ever appears to be much more complicated than the earlier one-dimensional 
problems. Figure (3) plots \5E/E \ for n > 2. 

For higher angular momentum states, there is a clear difference between 
the n = I — 1 cases and For the former case we find 5E oc —1/a 2 , a 

behaviour typical of nodeless states, while for latter case we find the expected 
5E oc —1/a trend. We explain the distinction between the two cases in 
Sect. (4). 

As for the dependence on the magnetic quantum number m we do have 
the expected invariance under m — > — m, but also find a mild dependence 
of the energy shift for for different m corresponding to the same n, I. For 
example, on a log-log plot of 5E vs. a for the n = 3, I — 2, m — 1, 2 states 
of the Hydrogen atom, the m = 2 line has slope —1.979(7) and intercept 
— 10.14(3) while the m = 1 line has slope —1.986(4) and intercept —9.61(2). 



4 General Analytical Investigation 

Unless otherwise stated, in this section we discuss the nonlinear equation in 
the presence of a general smooth external potential V(x) which for conve- 
nience we choose to be parity even, V(x) = V(—x), so that the unperturbed 
states are parity eigenstates. Since most studied potentials are parity even, 
that restriction is not unreasonable. However we emphasize that the key 
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features of our result, such as eg. (1441) below, follow from the structure of the 
nonlinearity ([2]): For example, without using parity eigenstates below, one 
still obtains similar results if instead of ([2]) one uses the L —>■ —L symmetrised 
version [5] of the nonlinearity. 

4.1 Nodeless States 

If p(x) does not vanish in the region of integration, such as the ground state 
of the SHO, one may use (TTOT) to conclude that SE ~ 0(L 2 ). Explicitly, we 
have for the n = SHO state, 



which for 77 = 1 /2 is in excellent agreement with the leading result extracted 
numerically in Section (3.2), SE = — 0.0156L 2 /a 2 . Equation fl25|) indicates 
a number of interesting features: It vanishes both as 77 — > 0, which is the 
formal linear limit of ([T]) and also as L/a — > which is the physical linear 
limit. SE also vanishes at rj = 1/3 and rj = 1 but it is apparent from ffTOl) 
that unlike the rj — > case the other two critical values are dependent on 



We also note that 5E in (1251) is positive for rj < 1/3 and negative for larger 
values. Such crossing behaviour will also be seen below for excited states but, 
more remarkably, at a universal (that is, V(x) independent) value of rj. 

The conclusion SE ~ O(L) 2 that we have drawn for nodeless states from 
({TO]) is for smooth one-dimensional potentials. For higher dimensions the 
conclusion is still true because of the form of (122]) but now one may encounter 
some nodes that are integrable, as in the hydrogen atom case to be discussed 
in Sect. (4.3) below. 

4.2 Excited States in One Dimension 

When the unperturbed wavefunction 4>(x) vanishes at a number of nodes the 
formal L-expansion of the quantum potential Q n i as in (flOl) breaks down and 
so one has to proceed differently. Now, from the definition of 5E in eq.(!7j) 
we have 



The second integral in (1261) is independent of L and it will cancel the L- 
independent piece of the first integral. So let us focus on the first integral 




(25) 



V(x). 




(26) 
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in (}2T)j) : Suppose first that p(x) has exactly one node at x = x\. Since 
there are two widely separated length scales, \L\ <C a, we may divide the 
integration region in the first term of (|26p into three parts, (—00, x\ — 2^)1 
[xi — x\ + ^y^] and (x\ + 00), where a is a positive constant to be 
fixed later. The absolute value \L\ used here allows negative L values in the 
following discussion. 

In the region including the node one may perform the Taylor expansion 
4>(x) ~ Ci(x — xi) and so p(x) ~ C 2 (x — £i) 2 . Thus 

tj node 

U L>i / 2 7 9 r i 

' dxx In 



4mL 2 ?7 4 JzBB. (1 - rj)x 2 + rj(x + r]L) 2 

(1 — rj)x 2 r)(x — r)L,y 



[1 — 7])x 2 + 7](x + r/L) 2 (1 — t])(x — r/L) 2 + r/x 2 



(27) 



ft C 2 |L| p/ 2 y 2 

— — / dyy 2 [ In — — + 1 

Ami]' 1 J -a/2 (1 — njy 2 + r)(y + r]) 2 



(1 - ??)y 2 + v(y + v) 2 (1 - v)(y - v) 2 + vv 2 

(28) 

Notice that the leading |L| 3 part of the integral (1271) comes already from the 
/ dxx 2 piece after the scaling x = \L\y, so subleading terms in the Taylor 
expansion of p n (x) ~ C 2 p (x — x p ) 2 + 0(x 3 ) contribute only at 0(L 2 ) to SE. 

For <p(x) having nodes at x — xi,x 2 , ....,x^ we may repeat the above 
procedure in the neighbourhood of each node as long as the nodes are widely 
separated. Then 



ti 2 \L 



N 



Ami] 

with 



SE nodes » ^4 J( V , a) J2 OL + 0(L 2 ) (29) 



np 
p=l 



J (77, a) = 



v)y 2 + v(y + v) 2 
(i-v)y 2 



+ 1 



viv-vf 



(1 - r])y 2 + r](y + r]) 2 (1 - r])(y - r]) 2 + -qy 2 



(30) 
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In (1251) n refers to the quantum number(s) of the unperturbed state and p 
labels a node. 

To fix the value of a we have to look at the nodeless regions of the the 
first integral in (1261) : The terms p(x ± rjL) may be safely expanded about 
L = to give for the integrand a series ~ L° + L 2 + there is no L l 
term in the series because of parity invariance (jSJ). But since the integration 
limits are dependent on a\L\ we have to be sure that the integral receives no 
enhancement of 1/\L\ factors from them and therefore we need to choose 

a= [L\ ) 

where C is a positive constant, so as to make the integration limits L- 
independent. Then, since SE = for L = 0, the 0(L°) piece from the 
nodeless regions of the first integral in (I26p must cancel the second integral 
in (126]) leaving a net contribution of order L 2 . 

Returning now to (129]) and using (13T|) . we deduce that for small values of 
our perturbative parameter \L\/a we need to expand J(rj, a) for a large. We 
find 



jfa a _> oo) _> _ _ rj rfl* ( 4?? - \y + O ^- = . (32) 
Define the a-independent piece of (1521) as 

so that one may finally write for (1261) 

*2|r| AT 

^=4^ J(7/) p E 1 C - +0(L2) - (34) 

A remarkable aspect of the formula (I34p is that the specific dependence 
on the external potential V(x) has been factorised: it is only in the C np 
coefficients. Since J(n) vanishes at 77 = r] c — 1/4, it means that the leading 
energy shifts vanish at a universal, V(x) independent, critical point! 

Given the intricate steps in above derivation, it is useful to perform some 
checks. From our numerical investigations we found that for the excited 
states of the infinite well, SE = for a value of rj between (0.24, 0.25), in close 
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agreement with the above prediction. We also confirmed numerically that 
the energy shifts for the excited states of the SHO also vanish at essentially 
(limited by our numerical precision) the same rj c , 0.24 — 0.25 as that for the 
infinite well. We remark that since the infinite well may be thought of as the 
7 — > oo limit of the potential V(x) = |x/a| 7 , so we are essentially checking 
( 13~4"|) at two limiting ends of a class of potentials. 

More complete checks of f[3~4"|) involve comparing it with the n-dependent 
expressions for the energy shifts found numerically in the previous section. 
These checks are done in the Appendices, again showing good agreement. 

Let us summarise some of the main features of (134"]) : Firstly, the expression 
clearly shows the non-analytic 0(|L|) trend confirmed numerically in the 
previous section. It also shows that the energy shift will be negative (positive) 
for large (small) rj values. 

4.3 Higher Dimensions 

For higher dimensions an explicit analysis similar to the preceeding subsec- 
tion is awkward because the nonlinearity is expressed in the preferred Carte- 
sian basis with broken rotational symmetry whereas most potentials, such as 
the hydrogen atom, have a symmetry and so are better expressed in other 
coordinate systems. Nevertheless, we can make some general statements. 

For nodeless states we have the analog of QUI) by expanding (I2"3"j) and so 
get SE ~ 0{L 2 ). 

For excited states the presence of nodes leads to singularities as before in 
the naive Taylor expansion. Arguments similar to above then imply that SE 
will be enhanced to 0(|L|) as each coordinate is treated separately in fl23|) . 
Thus we expect again the energy shifts to be positive for small rj and negative 
for larger rj, vanishing at some intermediate value. The numerical results of 
Sect. (3.3) for the angular momentum states / ^ n — 1 are in agreement 
with these general expectations although we have not checked the expected 
variation with rj. 

A very interesting situation arises for the I — n — 1 states of the hydrogen 
atom for which the radial wavefunction vanishes only at the origin, 

Vv=n-i ~ r^lUM) • (35) 

Although the corresponding probability density p(f) has a node at the origin, 
the radial integral in the 0(L 2 ) contribution SE ~ / dQ J °° drp(r) (...), the 
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three dimensional analog of (1101) . is nonsingular, as we see by power counting, 
if 2(n— 1)+3 > 4, that is, for n > 2. This explains the "anomalous" , 5E ~ L 2 , 
behaviour of such excited states observed in Sect. (3.3). 



5 Exactly solved models 

In using perturbation theory we have assumed, as is usually done in physics, 
that the quantity of interest will deform continuously as the perturbation 
is turned on. Here we briefly discuss two nonlinear Schrodinger equations 
for which exact solutions are available so that one can test perturbation 
theory. In addition, the models will be used to further highlight some of 
the distinctive features we have observed for the nonlinear equation (pQ). For 
simplicity we consider only the one dimensional case here. 



5.1 Gross-Pitaevskii (GP) Equation 

This classic [10] equation is used as an effective theory in studies of condensed 
matter. It corresponds to using F(p(x, t)) = gp(x, t) in (JT]). To leading order, 
one has 

SE = g J p 2 dx (36) 

so that energy shifts are always positive or negative depending on the sign 
of the coupling g. Explicitly for the infinite well one has 

a constant shift independent of n, as obtained earlier in Ref.[i~3] which also 
showed that this perturbative result was the appropriate limit of the exact 
solution of this equation with the infinite well potential. 

For the SHO potential we are not aware of any exact solutions for the 
GP equation but (I3"6T) gives the leading order correction, 

5E^^=n-°- 3 \ (38) 
ay 2tt 



for n > 1, showing that it decreases with n. We obtained fl38l) through a 
numerical best-fit to an assumed power law. 

The constant or decreasing dependence of the energy shifts on n, respec- 
tively for the above two potentials in the GP equation, should be contrasted 
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with the results for the information-theoretic nonlinearity (J2J) which showed 
an increasing dependence on n. As we saw in Sect. (4) that increasing depen- 
dence on n was due to the prominent role played by nodes which by contrast 
are completely irrelevant in fl36l) . 



5.2 A Pseudo-nonlinear model 

Starting from the usual linear Schrodinger equation 

h 2 d 2 ^ 

we can re-arrange the kinetic term by an amount e to obtain 

3^ , , Ti 2 d 2 ^ , , T e h 2 fd 2 ^\ . , 

ih— = -(l-e) — — r + M- T — hr 40 

which corresponds to an equivalent nonlinear Schrodinger equation with mass 
m/(l — e) in the linear part and a perturbed nonlinearity 

e ft 2 /3 2 *\ 

F(p) s " (41) 

Thus in this case the exact and unperturbed solutions just correspond to a 
mass renormalisation. For stationary states first order perturbation theory 
gives 

SB = /<i*J^ (42) 

2m J y/p ox 2 

2m J \ ox 

so that again the energy shifts are simply correlated in sign with the sign of 
e. 

As the wavefunctions for the linear Schrodinger equation with an infinite 
well potential are independent of the mass and so also of e, the first order 
correction using fT43|) gives an exactly e contribution in this case. When that 
is added to the unperturbed energies which are proportional to 1 — e, one 
gets the full answer, that is, first order perturbation theory for this problem 
is all there is. 



At <*m (43) 
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For the SHO and other problems the first order correction will generally, 
by construction, lead to final results correct up to errors of 0(e) 2 . 

Thus in this nonlinear model the perturbative corrections to the energy 
always have the same n dependence as that for the unperturbed energies. 
Therefore one may interpret the analogous results for Eq.(j2J) as due in some 
rough sense to higher derivative terms, higher than the second-order kinetic 
energy terms like (|41|) . This is indeed what is implied by a formal expansion 
of (|2J) but as we discussed in Sections (2,4), that formal expansion is in general 
singular and the actual result depends acutely on whether the unperturbed 
states do or do not have nodes. 

6 Conclusion 

Our main result is 

*2|r| N 

5E = -UL y/ v (l- V ) (1 - 417) £C 2 p + 0{L/af , (44) 

which gives the leading correction, due to the nonlinearity (j2J), to the energy 
eigenvalues of the usual 1 + 1 dimensional linear Schrodinger equation for 
cases where the unperturbed states have nodes. The correction is propotional 
to \L\/a, hence it is non-anayltic and an enhancement over the correction for 
states without nodes for which 5E oc L 2 . The dependence of 8E on the 
external potential is only through the C np coefficients defined in Sect. (4.2). 

From (j44p we see that independent of the external potential, V(x), the 
leading energy correction to states with nodes vanishes at rj = 1/4. The 
existence of such a universal critical point is quite unexpected as neither the 
equations of motion nor the full expression (j7J indicate such a special point. 

As eq. (!44p shows, 5E < for rj > 1/4, being positive for smaller rj. Since 
7] is a free parameter in the nonlinear equation (JTJ), it means that there is a 
qualitative difference in the properties of that equation for rj small or large. It 
is also interesting to note from (I44p that the expression is real precisely in the 
range < rj < 1, which is exactly the explicit condition on r\ we started with. 
Since for r\ — > one formally has the linear theory, the square-root factor 
again emphasizes, in addition to the \L\ term, the non-analytic character of 

dm). 

For the usual linear Schrodinger equation in one space dimension, states 
with nodes are the excited states of a system although in some cases, such 
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as that for the infinite well, the ground state also has nodes. For states 
without nodes, which are typically ground states, such as for the SHO, the 
leading energy corrections are of order L 2 and given by a simple expansion 
of ([7]). Thus at i] = 1/4 all the states of system, with or without nodes, have 
SE oc L 2 . 

For higher dimensions, the qualitative properties are similar to the one 
dimensional case. That is, nodeless states get SE oc L 2 while states with 
nodes in general have SE oc L. We saw an exception in the hydrogen atom 
example where some excited states with nodes had SE oc L 2 because the 
potential singularities were integrable. 

Let us now discuss the validity of perturbation theory for the infinite well 
and SHO where SE increases with n, the principal quantum number, faster 
than the unperturbed states. For example, in the infinte well case we found 
SE oc \L\n 3 /a so that even if \L\/a << 1, at large n the correction SE would 
no longer be small. This indicates a breakdown of first order perturbation 
theory for large n states, requiring one to go to higher orders. Presumably, 
if L/a is small, the net perturbative correction should be small for all n, so 
one expects the higher order corrections to sum to a reasonable expression. 
A simple Pade resummation suggests SE oc n 2 /(l + bn) for the infinite well 
at r] = 1/2, where b > is some constant. 

Finally, we discuss some physical implications of our result if the non- 
linearity (121) is a fundamental or effective representation of potential new 
physics at short distances as suggested in [5]. For a particle in a large box, 
we may use the infinite well result, generalised via (1221) to three independent 
dimensions, to see that for 7] > 1/4 high energy states have their energies low- 
ered, that is the nonlinearity acts to moderate high energy divergences. One 
reaches the same conclusion from the SHO results if one thinks of ordinary 
free quantum field theory modes as SHO states. 

Thus the nonlinearity d2J), applied here heuristically to field theory, sug- 
gests that the usual high energy divergences of quantum field theory might 
be moderated, if not absolutely eliminated. Now in [31 [5] it was suggested 
that the nonlinearity fl2J) might be linked to gravity simply by requiring L 
to be a universal length scale. Taken together, this then suggests that grav- 
ity might moderate ultraviolet divergences of field theory. Interestingly, the 
suggestion that gravity might regulate ultraviolet divergences has been made 
several times in the past through different reasoning within the context of 
usual linear quantum theory, see for example [14j and references therein. 

However the above moderation works only for rj > 1/4 where we have 
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SE < for excited states. What if in reality one has rj < 1/4? Then SE > 
and this means that we are quite possibly under-estimating the amount of 
energy in quantum systems. One wonders if this might be relevant for the 
dark energy/matter problem in cosmology. 

So it appears that knowing the physically relevant value of rj is quite im- 
portant for potential phenomenological applications of the nonlinear equa- 
tion. Perhaps rj could be fixed theoretically through a renormalisation group 
study of a discretised version of the nonlinear equation ([2]). In this regard, 
the naturally induced discretisation noted in [9] might be useful. 
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7 Appendix A: Infinite Well Revisited 



For the infinite well we may evaluate 5E (134"I) explicitly since the C n p 
these case are easily obtained from the wavefunctions (jTTj) . 



for 



a 



np 



2 nn , 

{-l) p , 0<p<n. 



a a 



(45) 



Thus the dimensionless energy shift is 



5E 



a 2 \L 



2vrV 



-n 



E 

P =i 

T] 4 



' 2 - 


f nir\ 2 


a 





(46) 
(47) 



The formula is valid also for the ground state, n — 1, because the correspond- 
ing wavefunction vanishes at the two end points, each of which contributes 
the equivalent of half of one regular node as can be seen by reviewing the 
derivation of (|29|) above. We therefore now have an understanding of the 
intriguing n 3 behaviour seen numerically in Sect. (3): each C 2 p contributes 
an identical n 2 piece to the sum over n terms. 
We find at r] = 1/2, 



5E n 



^16n 3 J(l/2) 
a 

\L\ , 
-1.05— n 3 , 
a 



(48) 
(49) 



in good agreement with our numerical study of the infinite well in Sect. (3) 
which indicated an average value of 1.03 for the numerical factor. 

For other potentials an explicit evaluation of the sum in (|29|) does not 
appear feasible as the coefficients in general are very complicated functions 
of n and p that are rarely known in closed form. However an asymptotic or 
numerical evaluation of C 2 p might be possible if an explicit dependence on 
n is required. We illustrate this for the SHO in the next Appendix. 
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8 Appendix B: Semi-analytical Analysis of 
SHO Energy Shifts 

Recall that our analytical estimates of the energy shifts for the excited states 
gave 

^=3 j( " ) g c - +o(i2) - 

For the SHO, we obtain C np from the wavefunctions 

^ n (z) = iV n (7ra 2 )- 1 / 4 F n (^)exp(-^/2), 



where z — x/a , N n = l/y/2 n n\, and H n (z)are the n-th order Hermite poly- 
nomials. Observe that the wavefunction vanishes only when the Hermite 
polynomial is zero. 

What is required is the Taylor expansion of the wavefunction about the 
nodes. At the nodes, 



H n ( Z ;) = for p=l,2,...ra, 



(50) 



where z™ refers to the p-th root of H n (z). Therefore near a root we have, to 
leading order, 



Reverting to x, one obtains 



2 I I a' l , n 



n {z p )exp 



(4) : '/2 



(51) 



Using the identity 
then gives 
Finally, 



H n+1 (x) + H' n (x) = 2xH n (x) 
H' n (z;) = -H n+1 (z%). 



(52) 
(53) 



19 



We evaluated the sum of these CyL through a numerical computation of the 
roots and sums of the Hermite polynomials. Since the 'a' behaviour is already 
explicit, we examined the n-dependence by calculating ^C\ v from n — 1 
to n — 23. Furthermore at 77 = 1/2: 

^ = T^(£^)^ = l/2) (54) 

= -0.27^n L40 . (55) 
a 

This result is in good agreement with the purely numerical one quoted in 
Section (3) (which was actually for relatively low values of a). 
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Figure Captions 



Figure 1 : Log-log plot of \SE\ vs. a for the ground state n — 1 of the 
infinite well. The line has slope 0.99985(2) and intercept 0.0448(1). In 
this and the other figures the 'logs' are natural logarithms. 



Figure 2 : Log-log plot of 



5E 



vs. n for a = 1000 of the SHO. The line 



has slope 1.413(8) and intercept -8.25(2). 

Figure 3 : Plot of \5E/E \ vs. n > 2 for / = states of Hydrogen atom. 
Curves for different values for a are shown. 
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